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Abstract 

We utilize the Shell Model Monte Carlo (SMMC) method to study the struc- 
ture of rare earth nuclei. This work demonstrates the first systematic "full 
oscillator shell plus intruder" calculations in such heavy nuclei. Exact solu- 
tions of a pairing plus quadrupole hamiltonian are compared with mean field 
and SPA approximations in several Dysprosium isotopes from A = 152 — 162, 
including the odd mass A = 153. Basic properties of these nuclei at vari- 
ous temperatures and spin are explored. These include energy, deformation, 
moments of inertia, pairing channel strengths, band crossing, and evolution 
of shell model occupation numbers. Exact level densities are also calculated 
and, in the case of 162 Dy, compared with experimental data. 

PACS numbers: 21.60. Cs, 21. 60.Ka,27.70+q,21.10.Ma 
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I. INTRODUCTION 



Our goal is to develop an improved microscopic understanding of the structure of rare 
earth nuclei; i.e., an understanding based on the behavior of individual nucleons in the 
nucleus. Toward that end we solve the shell model systematically in a full oscillator shell 
basis with intruders for the first time in rare earth nuclei using the Monte Carlo (SMMC) 
technique; calculations using other methods have been restricted to a severely truncated 
model space. SMMC allows us to trace structural rearrangements within nuclei induced by 
changes in temperature and spin, so that we may obtain a clearer microscopic picture of 
general structural features in this region of the periodic table. 

We assume an effective two-body nucleon-nucleon interaction and perform a Hubbard- 
Stratonovich transformation to obtain a path integral representation for the partition func- 
tion, which is then evaluated by Monte Carlo methods (see Section ||) to produce an exact 
shell model solution within statistical errors; this substantially enhances the predictive power 
of the nuclear shell model for some observables. Indeed, direct diagonalizations of the shell 
model Hamiltonian in a full basis have been limited to A~ 50, while we present calculations 
for A~ 150. 

We examine how the phenomenologically motivated "pairing plus quadrupole" interac- 
tion compares in exact shell model solutions vs. the mean field treatment. We also examine 
how the shell model solutions compare with experimental data; static path approximation 
(SPA) calculations are also shown. There have been efforts recently by others to use SPA 
calculations, since it is simpler and faster (see as examples). However, the SPA results 
are not consistently good. In particular, it is useful to know not only if phenomenological 
pairing plus quadrupole type interactions can be used in exact solutions for large model 
spaces, but also if the parameters require significant renormalization because this affects the 
accuracy of the SPA. 

We study a range of Dysprosium isotopes (Z=66, 86 < < 96), which exhibit a rich 
spectrum of the behaviors such as shape transitions, level crossings, and pair transfer that 
have been observed in the rare earths. These results should therefore apply quite generally 
in the rare earth region, although the immediate work focuses on Dysprosium. We have 
selected this element since the half-filled proton shell makes the model spaces particularly 
large. 

A previous paper discussed SMMC for the test case 170 Dy, which does not exist as a stable 
nucleus 0. The work presented here is much more systematic and thorough. Algorithm 
improvements subsequent to H have increased the computational execution by a factor of 
10 or more and have allowed us to calculate the rare earths at lower temperatures, nuclear 
shapes are calculated using the correct calculated quadrupole variance (not just a constant), 
and pairing operators not used in || are calculated. 

II. THEORETICAL BACKGROUND 

Shell model diagonalization is still limited to A ~ 50 in the Oflp shell [[|. In contrast, 
SMMC determines thermal observables, but explicit wave functions are never constructed; 
this is the key to how the predictive power of the shell model is extended so tremendously. 
The method is far less demanding on machine storage and there is no need to perform 
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manipulations with the exponentially increasing numbers of variables that are encountered 
in direct diagonalization. SMMC storage scales like N%N t , where N s is the number of single 
particle shell model states and N t is the number of time slices (see below). 

No known discrepancies exist between SMMC and direct diagonalization in cases where 
the comparison has been possible. This includes odd mass nuclei computed for appropriate 
temperatures. Realistic fp and sd-shell solutions using modified KB3 and Brown- Wildenthal 
interactions, respectively, agree with experiments []5| . These results give us a high degree of 
confidence in the SMMC technique. 

As with any shell model, an effective nucleon-nucleon interaction must be specified. 
We use the well-known pairing plus quadrupole interaction as formulated by Kumar and 
Baranger ||. The Hamiltonian is 

H = H sp - G p P}P p - G n PlP n - |Q • Q (1) 
with Q = Q p + Q n - The pairing and quadrupole operators are defined as 

jm 

Q-Q = J2 (3) 

i,j,k,l 

where = r 2 Y2^{9 , <ft) as usual. The single particle energies are also taken from Kumar 
and Baranger ||. 

Effective charges are incorporated to account for core polarization to fit measured electric 
quadrupole transition strengths. The electric quadrupole operator, with effective charges e p 
and e n , is 

Q & P Q P ~t~ &nQn (4) 



A. Method and sign problem 

Detailed procedures for the SMMC are explained fully in H and references therein. We 
provide no further explanation, except as regards the "sign problem." 

We define $ =TrC/ (T /|TrC/ er | as the sign for a given Monte Carlo sample, where U a is 
defined as 

U a = U Nt U Nt ^ 1 ...U 1 (5) 
Un = e~ A ^ (6) 

(7) 

h<j is the one body hamiltonian for the auxiliary field configuration a and (3 is the inverse 
temperature. In SMMC, the partition function path integral is divided into Nt time steps of 
size A/? so that (3 = N t A(3. Hence, the complete evolution operator is expressed as a product 
of operators in each time step. In these studies, we use the canonical (number projection) 
formalism to evaluate the trace. 
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If $ is not equal to one, numerical instabilities can arise. This has become widely known 
as the Monte Carlo "sign problem." The simple phenomenological pairing plus quadrupole 
interaction (without added pn pairing) does not have an inherent sign problem. However, 
sign problems can arise even with this simple interaction if time reversal symmetry is broken, 
as when odd masses are studied or the system is cranked by adding a term —ujJ z to H. In 
these studies, the sign violation turned out to be minor for odd mass ground states and 
canonical ensemble cranking was limited to uj < 0.3 MeV. Experimentally, these nuclei are 
observed to uj m 0.6 MeV 0. 



B. Shapes and moments of inertia 

The quadrupole expectation values vanish under rotational symmetry. However, 
for a given Monte Carlo sample will have some finite, nozero value. We calculate Qij = 
3xiXj — 5ijr 2 for each sample and relate its eigenvalues to the quadrupole (3 and 7 deformation 
parameters as done previously with SMMC f||3|]. The intrinsic frame for each sample with 
a field configuration a has nonzero components Q' and Q' 2 = Q'_ 2 as 

3 / 4 71 

(Q' 2 ). = |^\/f (r 2 ) CT ^sin 7 . (8) 
In terms of eigenvalues Q' u , Q' 22 , and Q' 33 of give 

(Q'n)a = \/f (V3«Q' 2 > CT + (Q'_ 2 ) a ) - V2(Q' } ( 



22). = (-V3«Q' 2 > CT + (Q'_ 2 ) a ) - V2(Q' ) a ) 



(Q 

We can also calculate the free energy F(/3, 7) to construct shape contour plots. This is done 
using 

FC9 )7 ) = -TlB^f (9) 
(j 6 sin 37 

P(/3, 7) is the shape distribution as a function of the deformation coordinates 7) and T 
is the temperature. Plots are truncated at small 7 for an obvious reason. All shape plots 
discussed in this paper are for the mass quadrupole and are not electric quadrupole. 

Moments of inertia are calculated in the cranked Hamiltonian (H — > H — uoJ z ) using 
X 2 = d(J z )/duj. We expect moments of inertia to initially increase as the nucleus is cranked 
and then to decrease when pairs are broken. 
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C. Pairing composition 



Pairing correlations in nuclei can be studied by calculating pairing strengths in different 
spins for protons and neutrons. For like particle pairs, define the pair creation operator as 

A ™ = x a iW-)' b (10) 

V 1 + a b 



where 

ft] a x &] b ]jM = (jam a j h m h \J M)a) a a) b (11) 

m a ,m b 



Note that we only consider like particle pairing, i.e., proton-proton and neutron-neutron 

(Eq. |l|). Now let a = (j a ,jb) and a' = (j c ,jd)- Using the pair creation operator, a matrix 
M^ a , can be constructed as 

= Y,{d JM (j a ,j b )Aj M (j c ,j d )) (12) 
hi 

from which we can then define a pairing strength P J as 

P J - = £ ( 13 ) 

a>a' 

The correlated pair strength, which is more useful, is obtained by subtracting uncorrelated 
mean field pairs from the total P Jw defined above. A Fermi gas has generally been used 
for the mean field with SMMC. Letting = (a\hk) and substituting nin 2 (5i3^24 — ^23^14) 
for {oia^a^a^) in Eq. [l^ yields the Fermi gas mean field pair strength P^ F . In this case, of 
course, we could use the SPA occupations as the "mean field" to subtract from the complete 
pairing plus quadrupole solutions. 

Even-even nuclei have correlated ground states, so we expect an excess of J = pairs 
beyond the mean field in even-even ground states. The hallmark for a pair condensate is 
the existence of one eigenvalue of M^ a , which is much greater than all the rest. 

The pair matrix can be diagonalized to find the eigenbosons B' aJM ^ as 

= ^i) a jMv{ab)A ] JMn (ab) , (14) 

ab 

where a = 1, 2, • ■ ■ labels the various bosons with the same angular momentum and parity. 
The ipajn are the eigenvectors of the diagonalization, i.e. the wavefunctions of the boson, 
and satisfy the relation 

ECj,V = ^- ( 15 ) 

jajb 

These eigenbosons satisfy 

^(BgJMTt-B'yJMn) — n aJ-K&a~i , (16) 
M 

where the positive eigenvalues n a j v are the number of J7i~-pairs of type a. 
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D. Backbending 



We can monitor the pair strength for neutrons coupled to J = 12 as a signature for the 
anticipated band crossing or backbending. The only orbital in our model space which can 
produce this coupling is the neutron z 13 / 2 level. We do not monitor backbending by mapping 
out the typical backbending plot of I vs. lo because the backbend in the plot requires a 
multivalued solution of I, whereas SMMC always produces a single- valued solution from the 
statistical ensemble. 



E. Level density in SMMC 

SMMC is an excellent way to calculate level densities. E{(3) = (H) is calculated for 
many values of f3 which then determine the partition function, Z, as 

\n[Z((3)/Z(0)] = - / VW); (17) 
Jo 

Z(0) is the total number of available states in the space. The level density is then computed 

as an inverse Laplace transform of Z. Here, the last step is performed with a saddle point 
approximation: 

S{E) = PE + \rZ{P) (18) 

p(£) = (2vr/r 2 Cr 1/2 exp(S), (19) 

where f3~ 2 C = —dE/d(3. SMMC has been used recently to calculate level densities in iron 
region nuclei |J. Here, we present the first exact level density calculation for the much 
heavier Dy. 

Nuclear level densities in the static path approximation has previously been investigated 



by Alhassid and Bush for a simple solvable Lipkin model ||10|| . The simple Lipkin Hamiltonian 
does not include pairing, however. These authors found SPA to be superior to the mean 
field approximation and the difference between the two depended on interaction strength. 



F. Recap on interaction and model space 

It is fortunate that the elementary pairing plus quadrupole interaction does not break the 
Monte Carlo sign so that no g-extrapolation is required. The calculation of level densities 
is also simplified, since accurate results require very low statistical uncertainties 0. 

We do not include isovector proton-neutron (pn) pairing in the interaction. This is 
a reasonable assumption since fp shell calculations with SMMC have shown clearly that 
isovector pn correlations diminish quickly as N exceeds Z and we are not at N~Z |TT 



Further, in our case the valence protons and neutrons occupy different oscillator shells. To 
be sure some pn correlations are included via the isoscalar quadrupole interaction (2Q P -Q n ), 
but any observable that depends strongly on pn correlations such as the Gamow- Teller 
strength will not be accurately determined with this interaction. 
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We choose one shell each for protons (sdg) and neutrons (pfh) with the opposite parity 
intruders hn/2 and 213/2, respectively. The space encompasses 32 proton levels and 44 



neutron levels. The oscillator length, b = yfi/mujQ, was taken to be 1.01A 1 / 6 fm. 

With this space, Dy has sixteen valence protons so that the proton shell is half-filled. 
The number of valence neutrons varies from four in 152 Dy to fourteen in 162 Dy. This model 



space is identical to the one used by Kisslinger and Sorensen |T2[, but is smaller than the 
two shell space Kumar and Baranger used. 



G. Static path approximation and mean field 

The static path approximation (SPA) is the one-time-slice limit of the partition function 
path integral and is obviously easily implemented in SMMC. The mean field approximation is 
the saddle point estimation of the path integral; neither includes imaginary time-dependent 
terms. The SPA differs from mean field since the integral over time-independent auxiliary 
fields in SPA is done exactly, so that contributions from multiple configurations (even with 
large fluctuations) are included rather than just the single, steepest-descent mean field. 

Static path calculations do not do as well with realistic interactions as with schematic 
interactions, e.g., pairing plus quadrupole. Also, the accuracy of the approximation varies 



among operators, as will be demonstrated in Section |HL 

Static path calculations were done in SMMC with continuous fields rather than the dis- 
crete field sampling that has typically been used with SMMC. However, continuous sampling 
is not used for other calculations due to the significantly increased number of thermalization 
cycles it would require. Continuous fields proved necessary for accurate results in SPA; cal- 
culated energies were higher for discretized fields than for continuous fields and the difference 
could not be attributed to statistical uncertainties. 



III. RESULTS 
A. SPA vs. full solution: static properties 

1. Energy, spin, and quadrupole moments 

Comparisons of SPA and full path solutions for the static observables energy (H), spin 
(J 2 ), and mass quadrupole moment {Q 2 ) are shown for the experimentally deformed, tran- 
sitional, and spherical isotopes in Figs. [l](a)-(i). A few of the full canonical calculations 
do not extend quite as low in temperature as the SPA results due to numerical instability 
developing from multiple matrix multiplications. Error bars in these plots are smaller than 
the dot sizes and are therefore not shown. Also, note that the quadrupole moments are 
expressed as (Q 2 )/b 4 , where b is the oscillator length. 

The SPA energy is greater than the exact energy, except at very high temperatures for 
all three of these nuclei; at the lowest temperatures, the difference is a few MeV. The origin 
of this discrepancy will be discussed below (Section [111 A 2\ ). The difference between SPA 
and exact canonical ground state energies is 2.39 ± 0.15 MeV for 152 Dy and 2.59 ± 0.18 
MeV for 156 Dy, so there is not a significant discrepancy between the lighter and heavier 
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isotopes. Thus, SPA does not predict absolute energies accurately, but works well for relative 
differences. The partition function integral in SMMC is always divided into time steps of 
fixed size A/3, which is fixed at 0.0625 MeV -1 . As the temperature, 1/(3, increases the 
number of time slices in the exact partition function expression decreases. Hence, it is not 
surprising that SPA is more accurate for higher temperatures. 

Looking at (J 2 ) shows that the SPA calculations only cool to about J = 8 — 10 for even 
the lowest calculated temperatures while the even-even ground states are of course J = 0. In 
these canonical SMMC calculations, J = is not exactly reached even in the full canonical 
calculations because the thermal ensemble always includes contributions from higher energy 
states. An estimate for the f3 required for good filtering to the even-even ground states is 
(3 = l/E(2~l), where E(2f) is the measured energy of the first 2 + state. This (3 value varies 
from (3 = 1.6 MeV -1 for 152 Dy to (3 = 11.5 MeV -1 for 162 Dy as E(2f) varies from 0.614 to 
0.087 MeV. 

The thermal spin expectation, (J 2 ), can, in principle, be compared against the experi- 
mental spectra. However, one never experimentally knows all the states in a nucleus so such 
a direct comparison is difficult except at low excitation energies which are dominated by 
the well known ground state band. The difference in (J 2 ) between SPA and full canonical 
solutions at the lowest temperatures is 58 for 152 Dy, 74 for 154 Dy, and 91 for 156 Dy. Hence, 
the deviation is worse with increasing deformation in these isotopes. 

SPA works very well for the quadrupole moments. This result is also very robust, i.e., 
strengthening or weakening the coupling x beyond the nominal Kumar-Baranger value does 
affect the agreement between SPA and exact results. 

More information about the static quadrupole moment appears in Figs. ^|(a)-(d), where 
the proton and neutron quadrupole components are shown separately. The neutron 
quadrupole moment increases quickly from 152 Dy to 154 Dy and increases much more slowly 
after the onset of deformation in 156 Dy. The proton quadrupole moment, meanwhile, remains 
approximately fixed for all the studied isotopes. (Q 2 ) increases rapidly as A increases from 
154 to 162. Relative to (Q 2 ) it is only 8.5% (Fig. 0(a)), but the total quadrupole moment 
at = 6 MeV" 1 is 46% larger for 162 Dy than for 156 Dy (1221 b 4 vs. 834 6 4 ). (Q 2 V )/(Q 2 ) is 
roughly 10% for all isotopes, where Q v is the isovector quadrupole operator (Q p — Q n ). 

Another comparison between SPA and full canonical calculations is how quickly the 
solutions cool for various observables. For example, in 152 Dy (H) appears to have stabilized 
by (3 = 6 MeV -1 in both SPA and the full solution. The spin (J 2 ) and quadrupole moment 
(Q 2 ) also appear to have minimized near f3 = 6 MeV -1 in the SPA and exact calculations. 
Similar results are evident in the other nuclei, except for (J 2 ) in 156 Dy which is clearly 
decreasing still in SPA for (3 = 8 MeV -1 , the largest value for which the calculation could 
be done. 

2. Pairing energy and gaps 

Some insight can be had by looking at the pairing energy and gaps. The quadrupole 
energy,— 0.5x(Q-Q) in SPA agrees well with the full canonical solution, but the total energies 
in SPA shown above have clear deviations from the exact canonical results. The difference 
in (H) is due to the pairing energy. The accuracy of the SPA in the pairing interaction 
depends strongly on the pairing strength and is naturally better for weaker pairing. 
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The pairing energies and BCS pair gaps for 152 ~ 156 Dy are shown in Figs. |3|(a)-(f). The 
latter were obtained from the pairing energies by A% cs /G = H pair . The disagreement 
between SPA and exact calculations looks worse for protons, but recall that there are 16 
valence protons and just 4-8 valence neutrons in these isotopes. 

B. B(E2) and effective charges 

The reduced electric quadrupole transition strength, B(E2), is computed from 

B(E2) = ((e p Q p + e n Q n ) 2 ) (20) 

where e p and e n are effective proton and neutron charges. We have taken e p = 1 + x, e n = x. 
Results are shown in Table [I] and Fig. |J where it has been assumed that the total B(E2) 
calculated in SMMC is the same as B(E2;2f — > Of). Effective charges in column 3 and 
column 4 are fitted to measured B(E2) values. Typical effective charges in rare earths are 
approximately e p ( n ) = 2(1), so these values are in a reasonable range. 

Table H and Fig. |4| also show what B(E2) strength would be obtained from SMMC 
quadrupole moments with Kumar-Baranger charges. This illustrates how a small change in 
the effective charge can produce a comparatively large change in the B(E2). For 162 Dy, the 
7% difference in e p when using Kumar-Baranger charges leads to a 23% change in B(E2). 

The collectivity of a nucleus, and thus the B(E2) value, varies with the energy E(2f). 
The effective proton charge, e p , is plotted against E{2f) in Fig. ||. For the SMMC results, the 
neutron effective charge needs to be zero for spherical 152 Dy to avoid severely overestimating 
the B{E2\ 2f — > Of) strength as it is calculated above. The fitted effective charge is e p = 1.75 
for the deformed isotopes and it is intermediate for 1 Dy. This is a reflection of the fact 
that exact solutions for the mean field interaction yield lighter dysprosiums which are too 
deformed (Section [III F|) ; the effective charge should be constant. 

It should be noted that Kumar and Baranger did not calculate B(E2) values for spherical 
nuclei in the same way as for deformed nuclei, i.e., they did not take effective charges for 
spherical nuclei as e p = 1 + 1.5Z/A, e n = 1.5Z/A. For spherical nuclei, they combined 
phonon and rotational model properties (see p. 552). They used E{2f) = [C jB) 1 ! 2 and 
B{E2,0 -> 2) oc Z 2 Rl{BC)- 1 / 2 from the phonon model with the relation Q s oc B{E2f/ 2 . 
These E(2f) predictions are not very good (Table |TJ). 

C. SPA vs. full solution: cranking 

D. Sign limits on cranking 

Recall from Section |II A| that cranking degrades the Monte Carlo sign from unity and 
that calculations become impractical when the sign drops below 0.5. This is illustrated for 
both full canonical and SPA results for 156 Dy in Figs, ^(a)-(b). Error bars are not displayed 
in these figures since they are small; the statistical error in $ for /3 — 6 MeV -1 and uj = 0.3 
MeV is 0.04. Canonical cranking in this nucleus is limited to small frequencies for (3 = 8 
MeV -1 , as evidenced by the quick drop in sign from uj = 0.05 MeV to uj = 0.1 MeV. The 
canonical cranking is fairly good for f3 < 6 MeV -1 , especially (3 < 4 MeV -1 . SPA cranking 
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predictably has better sign properties. In SPA, 156 Dy can be cranked well out to (3 = 8 
MeV -1 , which is the approximate limit of temperature that can be reached in this nucleus 
without matrix stabilization. 

1. Energy, spin, and quadrupole moments again 

Energy, spin, and quadrupole moments at various temperatures are compared at different 
cranking frequencies in exact canonical and SPA. Calculations for 156 Dy appear in Figs. 0(a)- 
(f). Energy results for SPA at different cranking frequencies mirror the full canonical result. 
In the range 5 < (3 < 8 MeV" 1 , or 0.2 < T < 0.125 MeV, uj = 0.1 MeV lies at small 
excitation energy e above the to = baseline, uj = 0.2 MeV lies roughly 3e above the 
baseline, and uj = 0.3 MeV is excited by approximately 6e above oj = 0. 

The spin results reveal that (J 2 ) in SPA is higher than the exact canonical result until 
T > 1 MeV when uj < 0.2 MeV. However, for u = 0.3 MeV, SPA agrees well with the 
exact solution. The exact solution is not shown for (3 = 8 MeV -1 at this frequency due 
to numerical difficulties. The SPA (J 2 ) for uj = 0.2 MeV is very flat across the computed 
temperature range, 0.125 < T < 1 MeV. As with the exact solution, (J 2 ) decreases with 
rising temperature for uj = 0.3 MeV. 

Quadrupole results with cranking are similar to the uj = results in that the quadrupole 
moment does not change when the temperature is decreased below 200 keV at any frequency 
studied here. Note that the quadrupole results are plotted vs. uj for various temperatures. 
At the lowest temperatures, the quadrupole moment begins to decrease after uj = 0.15 MeV 
in the canonical case. It decreases after uo — 0.1 MeV in the SPA, however, uj = 0.15 
MeV is not computed there so it is difficult to say if the quadrupole moment is declining at 
frequency 0.1 or 0.15 MeV in SPA. SPA agrees very well with the exact solution for (Q 2 ) at 
all temperatures and cranking frequencies computed. 

2. Moments of inertia 

The J z variation with cranking frequency determines the moment of inertia. Results for 
156 Dy are displayed in Figs. |8](a) — (b) for both canonical and SPA cases. 

The moment of inertia, X 2 , for 156 Dy is 44.6 ± 7 h 2 /MeV in the exact canonical ensemble 
at J « 4 with (3 = 8 MeV" 1 and is 73.0 ±2 h 2 /MeV at same uj in SPA. At this temperature, 
T = 0.125 MeV, the SPA moment of inertia is 64% larger than the exact result. The 
experimental moment of inertia I2 = 40fr /MeV at J = 4, which matches the SMMC 
canonical result. Also, the rigid body moment of inertia for 156 Dy with (/?, 7) = (0.24,0) is 
73h 2 /MeV. This coincides with the SPA moment of inertia. 

E. Band crossing 

The pairing strength for J = 12 pairs in 156 Dy, which can be produced only from i 13 / 2 
neutron pairs, is shown in Fig. |8|(c)-(d). The strength Pj=\2 begins to increase quickly in 
the canonical case for (/3,u) = (8,0.15), which corresponds to J = 16 ± 1. In SPA, Pj=\2 
increases sharply beyond ((3,uj) = (8,0.1), which corresponds to J = 14. 
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The occupations of both the proton hn/ 2 and neutron z 13 / 2 intruder orbitals for 156 Dy 
are given in Figs. 0(a)- (b). The proton intruder occupation is comparatively stable over 
this same spin range at each temperature. However, it is clear that the i\3/2 occupation is 
increasing with spin, particularly for lower temperatures, as expected. The h u / 2 occupation 
number decreases slightly with temperature for all frequencies computed. Occupation shifts 
slightly to nd 5/2 . For = 6 MeV" 1 , the maximum spin J z corresponds to J m 32 and for 
(3 = 1 MeV -1 , the maximum spin is J ~ 28. Unfortunately for (3 = 6 MeV -1 , or T = 167 
keV, the Monte Carlo sign is reduced to 0.4 at the maximum spin shown (recall Fig. |]). At 
(3 = 1 MeV -1 , however, the sign is still very stable at 0.96. 



F. Shape vs. temperature and spin 

Nuclear shapes have also been computed to clarify how the shape varies with tempera- 
ture and spin. Temperatures and frequencies for these calculations are given in the figure 
captions. In all shape graphs, the /3-axis is radial and the other axis is the 7-axis. Results 
for 152 Dy at temperatures from T = 0.25 MeV to T = 2 MeV are shown in Fig. ITUl The 
nucleus becomes increasingly spherical for rising temperature. 

Shapes for 154 Dy and 156 Dy are shown in Figs. [TT] - P~5] . These were all produced from 
the exact canonical ensemble except Fig. |12|(f), which was produced with the static path 
approximation. Cranked contour plots, such as Fig. |l~3"l(d) for 156 Dy at (/?, uf) = (8,0.1), 
show the nuclei becoming increasingly gamma-soft with increasing spin. This is also true in 



SPA (Fig. |12|(f)), which was utilized in this case since the Monte Carlo sign for the exact 
calculation becomes too small to obtain useful results. There is no sign of oblate shape at this 
spin in 154 Dy, as predicted by Cranmer-Gordon, et al., using a Nilsson-Strutinsky cranking 



model ||13|| . However, the SMMC ground state deformation in Dy with these parameters 



is clearly too large. Ma, et al. claimed evidence for a return to some collectivity in 



154 Dy from spin 36 + -48 + . The shape plot for 154 Dy at J ~ 50 (Fig. |i~2"l(f)) appears soft. The 
SMMC B(E2) = 5 W.U. at this spin using the fitted effective charge. 

Note that with increasing A in these isotopes, the ground state deformation is roughly 
constant and the depth of the minimum increases. In fact for 156 Dy, the depth of the well 
is roughly the same as the fission barrier (~ 40 MeV) |i~5| . However, it can be questioned 
whether our shape formalism (Eq. |]) should be applied at such low temperatures (T < 0.25 
MeV) since in the limit T — > the free energy becomes everywhere degenerate and zero. 
The very low temperature results for 156 Dy and 162 Dy (not shown) do not coincide with the 
mean field. 

Previous publication of shape plots from SMMC results in gamma-soft nuclei using a 
pairing plus quadrupole Hamiltonian with quadrupole pairing Jl6| did not exhibit this. 
However, those nuclei are only weakly deformed ((3 ~ 0.05 — 0.15) while the Dysprosiums 
with A>154 are well-deformed ((3 ~ 0.3). For the Dysprosium shapes in this paper, the 
unexpected depth of the potential well is only evident for well-deformed cases. 

It is apparent from the above that the dysprosiums are all deformed in their ground states 
in the exact model calculation. However, 152 Dy is known experimentally to be spherical. 
Kumar and Baranger did not calculate 152 Dy, though they did calculate some other spherical 
isotopes. A shape plot has also been constructed from SMMC results in 140 Ba for comparison 
with Kumar and Baranger (Fig. WM. The SMMC result for 140 Ba agrees with Kumar- 
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Baranger; both calculations indeed show a spherical nucleus. This isotope has Z = 56 and 
N = 84. For the shell model space used, this becomes six valence protons and two valence 
neutrons so that, unlike dysprosium, the proton shell is less than half filled. The SMMC 
B(E2; 2 — > 0) is 5 W.U. using the effective charge fitted for the dysprosiums. Reducing 
the quadrupole coupling to half its mean field strength still yielded deformation (3 ~ 0.3 for 
156 Dy with a deep potential. However, reducing x to half the mean field strength returns 
152 Dy to a spherical distribution which fits the measured B(E2) strength with effective 
charges (e p ,e n ) = (1,0). 

The equilibrium shape was also calculated in 144 Ba for inverse temperature (3 — 4 MeV -1 . 
144 Ba has E(2f) = 199 keV |23[ and deformation (3 = 0.19 ||24|| . This nucleus proved to be 
extremely deformed {(3 ~ 0.45) in SMMC with the Kumar-Baranger interaction parameters, 
but with a deformation well not nearly so deep as for the A > 154 Dysprosiums (Fig. [15]). In 
this case, the potential was only about 2.5 MeV deep. Kumar and Baranger did not calculate 
this isotope, so direct comparison with them is not possible in this case. Kumar and Baranger 



also made no claims their model is valid for nuclei at such extreme deformation 22 



G. Odd A 

As mentioned previously, the odd nucleon in an odd mass nucleus violates T-reversal 
symmetry and can break the Monte Carlo sign, even with an interaction free from repulsive 
contributions. Results from 153 Dy are shown below in Figs. |16|(a)-(e). In this case, for our 
simple Hamiltonian the Monte Carlo sign behaves well and remains at 0.82 for canonical 
(3 = 10. Also, some of this reduction in sign may in fact be due to limits of numerical 
accuracy in the machine. 



1. Static observables 

Results for static observables are quantitatively similar to the even-even results. The 
energy difference between the full canonical and SPA calculations at T = 100 keV is 1.97±0.4 
MeV (Fig. ^(a)), which is a little less than the 2.49 ± 0.16 MeV canonical-SPA energy 
difference in neighboring 152 Dy and the 2. 15 ±0.06 difference in 154 Dy. This difference is due 
to different pairing energies in these odd-even and even-even isotopes. The discrepancy in 
(J 2 ) between SPA and full solutions is AJ 2 = 66 or AJ pa 8 (Fig. |16|(b)). The ground state 



spin for 153 Dy is (7/2)~ so that (J 2 ) = 15.75 and the first excited state is (3/2)~ at excitation 
E = 109 keV. Thus, the estimated f3 needed for filtering the ground state is reachable and 
(J 2 ) = 16 ± 4 in the SMMC canonical ensemble agrees very well with experiment. Again, 
the SPA quadrupole moment in 153 Dy is in excellent agreement with the full canonical 
calculation (Fig. [I6|(c)). The total pairing energy and BCS gaps (Fig. 0(d)- (e)) are similar 
to results in 152 Dy and 154 Dy (Fig. 0(a)- (e)). 

Occupation numbers for protons and neutrons in the canonical ensemble for 153 Dy appear 
in Fig. |TT|(a)-(b). As the temperature increases to T = 2 MeV, the proton occupation shifts 
only slightly to the highest orbitals. For neutrons, however, there is a clear rise in the i 13 / 2 
occupation. The occupation numbers are also compared for full canonical vs. SPA in 153 Dy 
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for /3 = 10 in Figs. |17|(c)-(d). These occupation numbers look very similar, though the 
agreement is slightly better for the protons. Pairing strengths are more revealing. 

The pairing strengths in both proton and neutron channels has also been computed 
(Fig. (Ul(a)-(b)). The sum of these eigenvalues, with no background subtraction, in p(n) J = 
+ channels is 2.30(0.39) for the exact canonical solution and 1.70(0.27) in SPA. These values 
are stronger in the full canonical than in the SPA, as would be expected from looking at 
Fig. 0(b). For the protons, the difference in the eigenvalue sum is mostly due to eigenvalue 
number 2, where the full canonical eigenvalue is more than twice the SPA result. These 
eigenvalues are otherwise distributed very similarly in the full and SPA results. A similar 
situation holds for the neutrons, where the first eigenvalue for the full canonical solution is 
more than double the SPA value. SPA simply does not produce the nuclear pair condensate 
revealed in the exact calculation. 



2. Cranking 

With cranking at (3 = 10 MeV" 1 , the sign for 153 Dy is 0.69 for u = 0.05 MeV, 0.52 for 
w = 0.1 MeV, and just 0.10 for uo = 0.2 MeV. Recall that the sign for (3 = 10 MeV" 1 in 
uncranked 153 Dy is 0.82. The moment of inertia, X 2 , is 49.6 ± 3 h 2 /MeV in the limit lo — > 
for the canonical calculation. 



H. Level density 

The level density results for 154 Dy are shown in Figs. |T^- |2D|. E(/3) points are calcu- 
lated at intervals of A/3 = 0.0625 to execute the saddle point inversion to the level density 
p(E) (Eq. |TJ|). The level density (Fig. 0(a)) is compared with a few parameterizations of 
backshifted Fermi gas formulas. The 154 Dy density is not directly compared with experi- 
mental data since no measurements are available. SMMC results are not as accurate for low 
temperatures or small excitation energies (E < 1 MeV) since numerical errors tend to be 
larger there. This is not a serious concern since the saddle point approximation itself is not 
really valid at the lowest energies anyway. For the lowest energies, density of states is best 
determined by simple state counting from known experimental levels. 

Three versions of Fermi gas density formulas are used. The first, labelled BBF with 



a = 19.25 and S = 1.0 in Fig. |19|(a), is the classic Bethe formula [T? 



The calculation for 154 Dy was done with a = A/8 = 19.25 MeV -1 and the energy is back- 
shifted as E — S for 5 = 1 MeV for an even-even nucleus. This formula happens to agree 
quite well with the SMMC prediction for the 154 Dy density for energies above 2 MeV. Notice 
that solutions to this formula will diverge as E — > 5 for positive 5, so the result is only shown 
down to an energy where the density formula yields a sensible result. 

Holmes, Woosley, Fowler, and Zimmerman (HWFZ) calculate backshifted level densities 
as flTHl: 
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P{E) = ^§(E - 5)-^ exp (2y/a(E - 5)) (22) 

The a parameters for HWFZ can depend on whether the nucleus is deformed or not. For 
154 Dy, 5 = 0.89 MeV and a = 22.28 MeV" 1 for spherical parameters and a = 20.05 for 
the deformed parameterization. The spherical HWFZ curve is always slightly low and the 
magnitude tails off to quickly below E = 2.5 MeV as compared with the SMMC result. 
HWFZ (spherical) is too small by a factor 2.5 at E = 10 MeV and too low by a factor of 
4 at E = 1.5 MeV. HWFZ (deformed), which has a smaller a parameter, is clearly a worse 
fit. It is worse than an order of magnitude smaller at E = 10 MeV and is six times smaller 
at E = 1.5 MeV. The typical BBF a parameter, A/8, is 19.25 MeV" 1 for A = 154. This is 
smaller than the HWFZ (deformed) density parameter and would make the fit even worse. 

Thielemann has modeled the parameters 5 and a slightly differently than HWFZ. In this 
paper, these are called HWFZ-T parameters. He has taken 8 as 

5 = A(Z, N) - 10/ A (23) 

with 

A(Z, N) = 12/ \TA even - even (24) 

= -12/VA odd -odd (25) 

= odd (26) 

He obtained the density parameter a from a fit to experimental densities at one neutron 
separation energy |H5],|2pJ. For 154 Dy, this gives 5 = 0.90 and a = 19.58. The HWFZ-T level 
density is somewhat lower than the calculated SMMC density in 154 Dy at all energies, but 
the slope agrees pretty well with the SMMC calculation. The HWFZ-T magnitude is lower 
by a factor 15 at E = 10 MeV for 154 Dy and a factor 6 for E « 1.5 MeV. 

Certainly for the case of 154 Dy, the most naive Bohr-Mottelson Fermi gas formula works 
much better than the more carefully developed parameterizations of HWFZ and Thielemann. 
This serves as an example of the utility of more realistic SMMC calculations to determine 
nuclear level densities. 

From the specific heat (Fig. |20|) and the known E vs. (3, the 154 Dy density calculation is 



expected to be valid up to 10 — 15 MeV excitation before finite model space effects set in. 
The specific heat will increase with increasing temperature. Eventually, however, the model 
space will become exhausted as the valence particles are all promoted as high in energy as 
possible within the finite space. The turnover point where C v stops decreasing is taken as 
the limit of validity for the calculation. An inert core is assumed here at all times. 

The SPA level density for 154 Dy has also been calculated and compared with SMMC 
(Fig. |l9|(a)). The SPA level density agrees well with SMMC for low excitation energies, but 
is consistently lower for energies above 4 MeV. Recall the SPA energy E vs. (3 (Fig. |l|(d)) 
agrees with the full SMMC at high temperatures, but never cools completely to the SMMC 
value for lower temperatures. Thus AE from T = to T = oo is smaller in SMMC and this 
difference of a couple MeV makes a perceptible difference in the level density. At E — 10 
MeV, the SPA density is smaller by a factor of 4. 

The heat capacity can be found in Fig. ^1[ This looks similar to the full SMMC calcula- 
tion except the magnitude of C v is smaller except for the lowest temperatures (highest (3). 
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The heat capacity has a sharp dropoff below E = 1 MeV for both SPA and full SMMC solu- 
tions. The heat curve implies that SPA should be valid for up to about 13 MeV excitation. 
However, the SPA density clearly diverges from SMMC well before this limit. 

Similar calculations are shown for 162 Dy in Fig. |l9"l(b). For the more deformed 162 Dy, 
the HWFZ-T formula works comparatively well as Fermi gas estimates go, but is still off 
by a factor 3 near E = 1 MeV and factor 1.4 near E = 10 MeV. HWFZ in 162 Dy is better 
than HWFZ-T at low energies, but is clearly worse at higher energies. It is within a factor 
2 of SMMC for E = 1 MeV and smaller than SMMC by a factor of 3 at E = 10 MeV. In 
contrast to 154 Dy, HWFZ fits very well for 162 Dy using a = A/8 = 20.25 MeV" 1 . The simple 
backshifted Bethe formula fails badly here, however, especially for higher energies. 

We determined that the level density calculation for this isotope is valid up to excitations 
of 15 — 20 MeV. This is slightly higher than the valid range for the density in 154 Dy. 

The comparison of SMMC density in 162 Dy with the Tveter, et al data is displayed 
in Fig. The experimental method of Tveter, et al. can reveal fine structure, but does 
not determine the absolute density magnitude. The SMMC calculation is scaled to facilitate 
comparison. In this case, the scale factor has been chosen to make the curves agree at 
lower excitation energies. From 1 — 3 MeV, the agreement is very good. From 3 — 5 MeV, 
the SMMC density increases more rapidly than the data. This deviation from the data 
cannot be accounted for by statistical errors in either the calculation or measurement. Near 
6 MeV, the measured density briefly flattens before increasing and this also appears in the 
calculation, but the measurement errors are larger at that point. 

The measured density includes all states included in the theoretical calculation plus some 
others, so that one would expect the measured density to be greater than or equal to the 
calculated density and never smaller. We may have instead chosen our constant to match the 
densities for moderate excitations and let the measured density be higher than the SMMC 
density for lower energies (1-3 MeV). 

Comparing structure between SMMC and data is difficult for the lowest energies due 
to statistical errors in the calculation and comparison at the upper range of the SMMC 
calculation, i.e., E w 15 MeV, is unfortunately impossible since the data only extend to 
about 8 MeV excitation energy. 

Level density information has also been calculated for the lighter nearer closed shell 
nucleus 140 Ba. This was done to investigate possible systematic differences in level densities. 
Its level density is shown in Fig. |23| and the specific heat in Fig. |24|. Unlike the dysprosiums, 
the calculated heat capacity curve in 140 Ba is very flat. 



IV. SUMMARY 



The work has systematically laid the groundwork for applying the shell model in rare 
earths. Previous applications have been plagued by severely truncated model spaces. An 
advantage of being able to explore exact shell model solutions in more expansive model spaces 
is to explain in fundamental ways behaviors such as band crossings and pair correlations 
which have been previously understood from phenomenological models. 

The static path approximation for this phenomenological pairing plus quadrupole model 
works well for calculating deformation and relative energy differences between ground states 
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of different isotopes regardless of deformation. Additionally, deformations are well deter- 
mined in SPA for quadrupole coupling strengths even a factor of three larger or smaller than 
the Kumar-Baranger mean field values. 

The SPA results for the pairing energy and pair gaps are not as good, however, and 
the discrepancy is worse for increasing pair strengths. SPA also overestimates the low spin 
moments of inertia. However, SPA does produce the z^'13/2 band crossing at the predicted 
spin for 156 Dy. SPA does not produce the ground state nuclear pair condensate and pair 
gap, hence the discrepancies in energy and moments of inertia. 

Deformations in the canonical ensemble with Kumar-Baranger parameters agree with 
both Kumar-Baranger and experimental results for isotopes tested that appear in their 
paper ||, but deformations calculated for some other nuclei do not. Also, the deformation 
wells in Dysprosiums with A > 156 are very deep at low temperatures, i.e., below T ~ 
0.25 MeV. There is some question whether the shape plots at such low temperatures are 
reliable. However, 152 Dy is clearly too deformed in its ground state and the calculated 
B(E2) strengths require a reduced fitted effective charge for A < 154, indicating that the 
light Dysprosiums are excessively deformed in their ground states for this Hamiltonian. 

ACKNOWLEDGMENTS 

J. White is very grateful to T. L. Khoo, K. Langanke, W. Nazarewicz, and P. Vogel for 
useful discussions. This work was supported in part by the National Science Foundation, 
Grants No. PHY-9722428, PHY-9420470, and PHY-9412818. This work was also supported 
in part through grant DE-FG02-96ER40963 from the U.S. Department of Energy. Oak Ridge 
National Laboratory (ORNL) is managed by Lockheed Martin Energy Research Corp. for 
the U.S. Department of Energy under contract number DE-AC05-96OR22464. We also 
acknowledge use of the CACR parallel computer system operated by Caltech and use of 
MHPCC SP2 systems. 



16 



REFERENCES 



[1] R. Rossignoli, Phys Rev C, 54, 3, 1996. 

[2] R. Rossignoli, N. Canosa, and J.L. Egido, Nucl Phys A, 607, 3, 1996. 
[3] D. J. Dean, S.E. Koonin, G.H. Lang, PB. Radha, and W.E. Ormand, Phys Lett B, 317, 
3, 1993. 

[4] E. Caurier, et al, LANL nuclear theory preprint #9809068, 1998. 

[5] D.J. Dean, S.E. Koonin, and K. Langanke. Phys. Reports 278, 1 (1997). 

[6] M. Baranger and K. Kumar, Nucl. Phys. A. A 110, 529 (1968). 

[7] Private communication from T.L. Khoo, 1998. 

[8] W.E. Ormand, D.J. Dean, C.W. Johnson, G.H. Lang and S.E. Koonin, Phys. Rev. C49, 
1422 (1994). 

[9] H. Nakada and Y. Alhassid, Phys Rev Lett, 79, 16, 1997. 
[10] Y. Alhassid and B. Bush, Nucl. Phys. A549, 43, (1992). 

[11] K. Langanke, D. J. Dean, S. E. Koonin, and P. B. Radha, Nucl Phys A, 613, 3, 1997. 

[12] L.S. Kisslinger and R.A. Sorensen, Rev. Mod. Phys. 35, 853 (1963). 

[13] H.W. Cranmer-Gordon, et al., Nucl. Phys. A465, 506 (1987). 

[14] W.C. Ma, et al, Phys Rev Lett 61, 1 (1988). 

[15] W. D. Myers and W. J. Swiatecki, Nucl. Phys. 81, 1 (1966). 

[16] Y. Alhassid, G.F. Bertsch, D.J. Dean, and S.E. Koonin, Phys. Rev. Lett. 77, 8, 1996. 
[17] H.A. Bethe, Rev. Mod. Phys. 9, 69, 1937. 

[18] Holmes, Woosley, Fowler, and Zimmerman, Atomic and Nuclear Data Tables 18, 305 
(1976). 

[19] G. Rohr, Z Phys A 318, 299 (1984). 

[20] J.J. Cowan, F.K. Thielemann, and J.W. Truran, Physics Reports 207, 1991. 
[21] T. Tveter, L. Bergholt, M. Guttormsen, E. Melby, and J. Rekstad, Phys Rev Lett 77, 
2404 (1996). 

[22] M. Baranger and K. Kumar, Nucl. Phys. A, A 110, 490 (1968). 

[23] Table of Isotopes, 8th ed. Richard B. Firestone and Virginia S. Shirley editors. John Wiley 

and Sons, Inc. New York, 1996. 
[24] Atomic and Nuclear Data Tables, 36(1), 7 (1987). 



17 



FIGURES 
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FIG. 1. (a)-(c) Energy, spin, and quadrupole moment for 156 Dy. (d)-(f) Energy, spin, and 
quadrupole moment for 154 Dy. (g)-(i) Energy, spin, and quadrupole moment in 152 Dy. 
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FIG. 2. SPA and exact canonical solutions for quadrupole observables in selected dysprosium 
isotopes A = 152 to A = 162. All results are for (3 = 6 or T = 0.167 MeV. {Q 2 V ) is denned as 
(Qp — Qv) 2 and b is the oscillator length. 
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FIG. 3. Pairing energy in 152 Dy (a), 15 Dy (b), and 156 Dy (c). BCS pair gaps are shown for 
152 Dy (d), 154 Dy (e), and 156 Dy (f). 
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FIG. 9. Occupation for proton hn/ 2 (a) and neutron i 13 / 2 (b) vs. spin. 
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(a) 152 Dy, T= 0.25 MeV (b) 152 Dy, T= 0.5 MeV 
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FIG. 10. (a) Shape at T=0.25 MeV for 152 Dy. Contour spacing is 0.79 MeV. (b) Shape at 
T=0.5 MeV for 152 Dy. Contour spacing is 0.64 MeV. (c) Shape at T=2 MeV for 152 Dy. Contour 
spacing is 2.9 MeV. Each plot is compiled from 2000 sample. 
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FIG. 11. Shape at T=0.133 MeV for 154 Dy. 4800 samples. Contour spacing is 0.92 MeV. 
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FIG. 12. (a) Shape at T=0.25 MeV and u = 0.05 (J » 6) for 154 Dy. 2560 samples. Contour 
spacing is 2.9 MeV. (b) Shape at T=0.25 MeV and u = 0.1 (J « 8) for 154 Dy. 2400 samples. 
Contour spacing 2.5 MeV. (c) Shape at T=0.25 MeV and to = 0.2 for 154 Dy. 2400 samples. 
Contour spacing 2.4 MeV. (d) Shape at T=0.5 MeV and uj = 0.05 (J w 10) for 154 Dy. 3840 
samples. Contour spacing 3.4 MeV. (e) Shape at T=0.5 MeV and ui = 0.2 (J rs 20) for 154 Dy. 
1920 samples. Contour spacing 4.2 MeV. (f) Shape at T=0.125 MeV and uj = 0.6 (J w 50) for 
154 Dy in SPA. Sign $ = 0.5. 2000 samples. Contour spacing 1.2 MeV. 
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FIG. 13. Shape at T=0.25 MeV for 156 Dy. 2000 samples. Contour spacing is 6.5 MeV. Shape 
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156 Dy. 2000 samples. Contour spacing is 2.0 MeV. Shape at T=0.125 MeV for 156 Dy with u = 0.1 
MeV' 1 . 4800 samples. Contour spacing is 3.0 MeV. 
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FIG. 14. Shape at T=0.25 MeV for 140 Ba. 4000 samples. Contour spacing 0.42 MeV. 
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FIG. 15. Shape at T = 0.25 MeV for 144 Ba. 2300 samples. Contour spacing is 0.7 MeV. 
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FIG. 16. Energy (a), spin (b), quadrupole moment (c), pairing energy (d), and BCS pair gap 
(e) in 153 Dy for Kumar-Baranger interaction strengths. Symbols for (e) match the symbols in (d). 
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Pairing matrix eigenvalues for neutrons in J 71 = + for 153 Dy at T = .1 MeV. 
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FIG. 19. (a) Level density for 154 Dy in SMMC shown with backshifted Fermi gas approxima- 
tions. BBF=standard backshifted Bethe formula. HFWZ=parameterization of Holmes, Woosley, 
Fowler, and Zimmerman, and HWFZ-T=HWFZ formula with Thielemann parameters. See text. 
HWFZ with a = 22.28 is for spherical 154 Dy parameters while a = 20.05 is for deformed 154 Dy 
parameters. SPA denotes the SPA result for the level density from SMMC. (b) SMMC result for 
162 Dy compared with the same backshifted Bethe formula approximations. 
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FIG. 21. SPA heat capacity in 154 Dy. 
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FIG. 23. SMMC density in 140 Ba. 
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FIG. 24. SMMC heat capacity in 140 Ba. 
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TABLES 



TABLE I. SMMC B(E2) vs. measured B(E2;2f -> Of) with specified effective charges. 
B{E2) in W.U. Errors are statistical Monte Carlo sampling errors. Column 6, B(E2) for 
SMMC/KB e p (n), is the B(E2) obtained from SMMC calculated quadrupole moments with Ku- 
mar-Baranger effective charges. 



A 


N 






B(E2) 


B(E2) 


B(E2) 










SMMC 


SMMC/KB e p(n) 


Expt 


152 


86 


1 





13 ±0.2 


N/A 


13 


154 


88 


1.5 


0.5 


97 ±0.2 


N/A 


97 


156 


90 


1.75 


0.75 


146 ± 0.6 


126.4 ±0.5 


146 


162 


96 


1.75 


0.75 


199 ±0.7 


150.1 ±0.5 


199 



TABLE II. Some spherical E(2f ) values from Kumar-Baranger 0. Energies in keV. 





^(2^ )ea;p 


E(2f) th 


138 Ba 


1438 


2767 


140 Ba 


602 


2006 


140 Ce 


1596 


2531 


142 Ce 


641 


1772 


144 Ce 


397 


1095 
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